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P^' Abstract 



A generalized quantum search algorithm, where phase inversions for the 
marked state and the prepared state are replaced by tt/2 phase rotations, is 
realized in a 2-qubit NMR heteronuclear system. The quantum algorithm 
searches a marked state with a smaller step compared to standard Grover al- 
gorithm. Phase matching requirement in quantum searching is demonstrated 
by comparing it with another generalized algorithm where the two phase ro- 
tations are tt/2 and 3n/2 respectively. Pulse sequences which include non 90 



degree pulses are given. 
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Grover's quantum search algorithm is one of the most important development in quantum 
computation ||1]]. It achieves quadratic speedup in searching a marked state in an unordered 
list over classical searching algorithms. It has many potential applications in various fields 
of interests, for instance in deciphering the DES encryption code and algorithms that 
need searching. Typical examples of the algorithms that need searching are the Simon 
problem P], the Hamilton's circuit problem f^, the hiddens shift problem [Q and quantum 
counting ^. There are several generalizations of the Grover algorithm. A modification 



of the algorithm can search for a "chain" of m marked items in 0{J N/m) iterations [0. 
In the standard Grover algorithm, the quantum database is built on an evenly distributed 
quantum superposition, and a generalization is made to allow the algorithm to work on a 
biased database where the amplitudes of the items in a database are not even p. Searching 
in an arbitrary entangled superposition is given in Ref. 

In some cases, one needs a quantum searching engine that searches an item with a smaller 
step. For instance, in the Simon algorithm 0] where ^-phase rotations rather than phase 
inversions are used and in the case where the number of marked states is more than N/4: ||10|| , 
standard Grover algorithm can not be used. In addition, for small A^, state of the quantum 
computer may not be exactly the marked state during the search process, and there is small 
probability that the algorithm may fail. In problems that certainty of success is vital, this 



should be avoided. A generalized quantum search algorithm |[TT|JT^ suits this purpose, where 
the searching step can be anything between that of Grover algorithm and zero. This is done 
by replacing the two phase inversions in Grover's algorithm with smaller phase rotations(0 
for phase rotation of the marked state, and 9 for phase rotation of the prepared state |0...0) 
). It has been found that with an evenly distributed database, arbitrary phase rotation 
is not applicable [|1T|, and only when the two phase rotations satisfy the phase matching 



requirement 6 = (p that an efficient quantum searching algorithm can be constructed [12 



An approximate expression 2sin| J j^ was given for the search step |T2|, where 6 is the 



rotation angle of the marked state [O. Exact expressions are given in a S'0(3)geometric 
picture [jl3 . 



NMR implementation of Grover's algorithm has been realized in 2-qubit and 3-qubit 
systems |I^ , p^ , p!7| . Quantum counting has also been realized in NMR system |18[. Exper- 
imental studies are important in demonstrating quantum algorithms, investigating effects 
of gate imperfection and decoherence, and in identifying problems in building a practical 
quantum computer. In this Letter, we report the experimental realization of this phase 
matching quantum search algorithm where the phase inversions are replaced by rotations 
through 7r/2 in a 2 qubit heteronuclear system using NMR technique. To demonstrate the 
effect of phase matching, another experiment where 9 = 7r/2, (p = 37r/2 is also performed. 
Different from standard Grover algorithm, the pulse sequences used here contain non 90 
degree pulses, and the delay pulses need not be multiples of jj. 

Grover algorithm consists of four steps in an iteration [^: l)a phase inversion of the 
marked state Ir = / — 2|r)(r|; 2) the Walsh-Hadamard transformation W; 3) a phase inver- 
sion of the prepared state |0), Iq = / — 2|0)(0|; and 4) the Walsh-Hadamard transformation. 
The operator for one Grover iteration is Q = —WloWIr- The steps 2-4 are combined to 
give the inversion about average D which has the following matrix 

f - i^i 
D,, = WIrW={ ^ ^' (1) 

In the generalized quantum search algorithm, the phase inversions are replaced by arbi- 
trary phase rotations. The corresponding operators, indicated by a "g" in the superscript, 
are 

Q9 = D9 J9 = ]Y J9 ly J9^ (2) 

where 
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When 6 = (p = IT the Grover algorithm is recovered. It is helpful to give the detailed 
expressions for a 2-qubit system. We assume that the marked state is r = 3. The operators 
are, 
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and 



Q9 = W I^W 
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For demonstration in this experiment, we choose 9 = cf) = 7r/2. In table |, we give 
the state vector for the quantum computer in each step during a quantum searching process 
where \c) = l/\NY,i^T K)- When phase matching requirement is not satisfied, the quantum 
algorithm will not work. To demonstrate this, we also give the state vectors of the quantum 
computer for the case with 6 = 7i/2 and = 37r/2. Both algorithms are performed for 
10 searching iterations. The probability for finding the marked state is the square of the 
coefficient of |r). These state vectors are converted to density matrices for comparisons with 
experiment. 

In the experiment, the working media is H2PO3. The 2 qubits are the nuclear spins of 
the H-atom and the P-atom. The observed J-coupling between ^H and ^^P is 647.451 Hz. 
The experiment was performed on a Bruker 500 ARX NMR spectrometer. The frequencies 
for ^i7(qubit A) and '^^P (qubit B) are 500MHz and 220MHz respectively. The Hamiltonian 
for this system can be modelled as a two-spin system with a weak coupling interaction, 



H = loa Iaz + ^B Ibz + '271 Jab Iaz Ibz, 



(6) 



where Iaz = \^za is the angular momentum operator in the e^ direction for spin A. uja 
and ujb describe the free precession frequencies of the nuclear spins A and B respectively. 
The magnetic field is in —el direction. The quantum gate operations needed in quantum 
computation can be constructed by a combination of some hard pulses and the delay pulses. 
Compared to the pulse sequences in the Grover algorithm for 2-qubit system, one needs only 
a small modification in the pulse sequences. First let's denote the following operators 

X^ = exp[-iei^], 

X^ = exp[i9Ix], 

Y^ = exp[-ieiy], 

Y^ = exp[ieiy], (7) 

which are radio frequency(rf) pulses for rotations about x-axis through 6, —6, rotations 
about y-axis through 9 and —9 respectively. In addition to these hard pulses, we also have 

r* = exp[-2TTi Jab Iza Izb t], 

which is a delay pulse where the system undergoes an evolution during period t in the doubly 
rotating frame. We denote | ]) = |0) and | |) = |1). 

We used temporal averaging |jl9l to produce the effective pure state |00). The pulse 
sequences for the temporal averaging and the Hadmard- Walsh transformation are standard 



14,15: 



Pq : I{none] 



PI : Fjr^X|Y2r^Xl, 



P2:17r^XiXir^Xi, 



W = {XjfYAHXlfYB^ . 



The pulse sequences for the generalized quantum search algorithm are obtained by modifying 
the pulse sequences used for Grover's algorithm in Ref. |T^: 






TT TT 



27r- 



D^ = YjXlYA^YiXlYB^T^y^. (8) 

The pulse sequence for a complete generalized search iteration Q^ = D^ P is 

Q^ = XA^17Xirir^i^Xir^^XirB^r^7^. (9) 

It is interesting to study the length of time in a given quantum search algorithm with 
arbitrary phases. A hard pulse takes microsecond to complete while a delay pulse takes 
about a millisecond to complete. An algorithm with large 6 or (p takes less time to complete. 
Therefore in practice, it is better to use large phase rotations to make the computation time 
short if one is given the freedom. It also worth pointing that in general, the hard pulses may 
not be the multiples of | as in other NMR quantum computations, and the delay pulses 
may also takes noninteger multiples of ^. This is different from the pulse sequences used 



in standard Grover algorithm, for instance in |14,17 



Two sets of experiments: phase matched searching with 6' = = | and phase mismatched 
searching with ^ = |,0 = ^, have been performed. State tomography is used to obtain 
the density matrices [^. Density matrices are experimentally constructed for all the 10 



iterations in the two sets of experiments. It took quite some time and labor to get them. In 
table H, we have given the relative error defined as 6p = *''|7^r^ • However, this error is 



not solely the "genuine" errors pT| , p2| from gate imperfection and decoherence that occurs 
during the quantum computation process. Part of the error is caused by doing the integration 
of areas of the spectrum by hand during the density matrix construction. It is interesting 
to note that the relative errors at later stages are sometimes even smaller than the early 
stages. For instance at step 7 in the phase-matched case, the relative error is only 15% , the 
smallest among the 10 steps. Similar result is observed in other NMR quantum computing 
experiment, for instance in |T^. This is perhaps because imperfect gate errors cancel out 



each other. For the economy of the paper length, we give only the density matrix for the 6th 
iteration where the success rate is the largest in Fig.l for phase matched searching. In Fig. 2 



the same is plotted for the phase mismatched searching. For clarity, we plot separately 
the real and imaginary parts. From these figures, it is seen that the agreement between 
theoretical prediction and experimental data is good. In particular, when phase matching is 
satisfied, the algorithm can find the marked state with high probability. However, when the 
phase matching requirement is seriously violated, the probability of finding the marked state 
is very low. Phase mismatching leads reduction in success rate, and it should be avoided in 
practice. However, Grover's algorithm has some intrinsic robustness against errors, a small 
amount of phase mismatching, like those errors from NMR pulse manipulations, will not 
cause a big loss in the probability for finding the marked state. 

In summary, a generalized quantum search algorithm has been demonstrated in a 2- 
qubit NMR system. Pulse sequences are given. Non 90 degree pulses are used and tested 
to give good performance. It also demonstrate that phase matching in quantum searching 
is important. 

This work is support in part by China National Science Foundation under Grant 
No. 60073009, the excellent university teacher's fund of China Education Ministry, Fok Ying- 
Tung education foundation. 



TABLES 





TABLE L Theoretical State Vector of the Register in Each Search Iteration 




phase matching 


phase mismatching 

= 7r/2,(/> = 37r/2 


stepl: 


\ipi) = 0.90|t) + 0.43|c) 


iVi) =0.25|r)+0.97|c) 


step2: 


|V'2) = 0.97|t) + 0.22|c) 


1^2) =0.62|r)+0.78|c) 


step3: 


1^3) = 0.65|t) + 0.76|c) 


IV3) = 0.06|r) + LOO|c) 


step4: 


IV'4) = 0.39|t) + 0.92|c) 


1^4) =0.59|r)+0.80|c) 


step5: 


|V'5) = 0.78|t)+0.62|c) 


IV5) =0.36|r)+0.93|c) 


step6: 


|V'6> = L00|t) + 0.01|c) 


iVe) =0.41|r)+0.91|c) 


step7: 


IV'7) = 0.80|t) + 0.60|c) 


1^7) =0.57|r)+0.82|c) 


step8: 


IV's) = 0.40|t) + 0.92|c) 


l^s) =0.13|r)+0.99|c) 


step9: 


|V'9) = 0.63|t) + 0.77|c) 


1^9) =0.63|r)+0.78|c) 


steplO: 


|V'io) = 0.97|r)+0.24|c) 


|V'io) = 0.19|r)+0.98|c) 



TABLE II. The errors of the experiments 



Stepl %18 %17 

Step2 %21 %28 

Steps %20 %27 

Step4 %21 %21 

Step5 %27 %20 

Step6 %20 %20 

Step? %15 %16 

Steps %24 %33 

Step9 %22 %30 

SteplO %22 %33 
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Figure captions 

Fig.l Comparisons of theoretical and experimental density matrices with 6 = (p = tt/2 for 

iteration 6. (a) and (b) are the theoretical real part and imaginary part of the density 

matrices respectively, whereas (c) and (d) are the corresponding experimental ones. The 

probability of finding the marked state is the sum of squares of the 11 component which are 

in the upper right corners of the figures. 

Fig.2 Same as Fig.2 but for 6 = 7r/2, = 37r/2. 
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